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We investigate the transmission and reflection survival probabilities for the chaotic stadium bil- 
liard with two holes placed asymmetrically. Classically, these distributions are shown to have al- 
gebraic or exponential decays depending on the choice of injecting hole and exact expressions are 
given for the first time and confirmed numerically. As there is no reported quantum theoretical or 
experimental analogue we propose a model for experimental observation of the asymmetric transport 
using semiconductor nano-structures and comment on the relevant quantum time-scales. 
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Billiards [Ij are systems in which a particle alternates 
between motion in a straight line and specular reflections 
from the walls of its container, while open billiards con- 
tain one or more holes through which particles may es- 
cape. Billiards demonstrate a broad variety of behaviours 
including regular, chaotic and mixed phase space dynam- 
ics, depending on the geometry, whilst allowing for math- 
ematical treatment of their properties. Billiard models 
have been increasingly important in both theoretical and 
experimental physics, for example as models in statistical 
mechanics such as the Boltzmann Hypothesis , number 
theory and the Riemann Hypothesis [3 , in room acous- 
tics [4 , atom optics, where ultracold atoms reflect from 
laser beams [5 , optics in dielectric micro-resonators [6] 
and in quantum chaos when solving the Helmholtz equa- 
tion with Dirichlet or Neumann boundary conditions [7] . 
Open billiards are also a useful model for understanding 
the close correspondence between classical and quantum 
mechanics 0. 

Quantum open billiards were experimentally realized 
first in flat microwave resonators in the early 90 's [9 , 10] 
and later in semiconductor nano-structures such as quan- 
tum dots [TTl |T2]. Experiments perturbing these sys- 
tems with small magnetic fields exhibit principal quan- 
tum interference effects like weak localization, Altshuler- 
Aronov-Spivak oscillations and conductance fluctuations, 
all of which semiclassical theory has arguably succeeded 
to explain using properties of the underlying classical 
dynamics [T2j [13]. Similarly, in microwave resonators, 
due to their clean, impurity-free geometry and the tun- 
able coupling strength to the various decay channels, pre- 
dicted phenomena like resonance trapping have been ex- 
perimentally observed [M]. 

Here we investigate the classical transport of a pop- 
ular example for the above and other experiments, the 
stadium billiard with two holes on its boundary (see Fig- 
ure [T]). Looking at the phase space of this open system, 
we find that the predominantly chaotic character of the 
corresponding closed system is non-trivially affected by 
the positioning of holes. In particular we find that the 
transmission and reflection probabilities, when particles 



are injected from one of the two holes, are qualitatively 
different at long times depending on the choice of the 
injecting hole therefore displaying time-dependent asym- 
metric transport. We give detailed analytical expressions 
for these distributions and confirm them numerically. Al- 
though some work has been done in this direction on the 
quantum level [15 , to the best of our knowledge, there 
has been no analogous analytic prediction or experimen- 
tal observation of such an effect. Hence we conclude with 
a discussion of a possible experimental model with re- 
gards to the relevant quantum time scales. 
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FIG. 1: Stadium billiard with two holes Hi and H2. The billiard 
map is parameterized using arc length < x < 4a + 27rr and 
velocity parallel to the boundary -usin^ with G (—^,^). The 
hole's coordinates on the straight segment are defined such that 
—a < h^ < hf < a. 



The transport problem is closely related to the escape 
problem for which we also make new observations. The 
uniform (Liouville) distribution projected onto the bil- 
liard boundary has the form {2\dQ\)~^dx dsin^ (where 
I^QI is the perimeter of the billiard while x and are 
defined in Figure [T]), and is the most natural choice for 
an initial distribution of particles. Given such a distri- 
bution, the probability P{t) that a particle survives {i.e. 
does not escape through k small holes Hi G dQ) in a 
strongly chaotic billiard up to time t decays exponen- 
tially ~ e~^^ at long times with the exponent to leading 
order given by [16 
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where (r) = is the mean free path for 2D billiards. 
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hi = \Hi\ is the length of each hole, \Q\ the area and v 
the speed of the particles. 

The stadium billiard is a chaotic system where the de- 
focusing mechanism guarantees a positive Lyapunov ex- 
ponent A (exponential separation rate of nearby trajecto- 
ries) almost everywhere [17j, the exception being a zero- 
measure family of marginally unstable periodic orbits be- 
tween the parallel straight segments called Bouncing Ball 
orbits. They have been shown to lead to an intermit- 
tent, quasi-regular behaviour which effectively causes the 
closed stadium to display some weaker chaotic properties 
such as an algebraic decay of correlations [TQ . Quantum 
mechanically they cause scarring [19j, the system is not 
quantum uniquely ergodic [20], an h dependent 'island of 
stability' appears to surround them [ST and deviations 
from random matrix theory (RMT) GOE predictions are 
observed (especially in the Aa-statistics) if not treated 
appropriately (see [22 , 23 ). 

A small hole of size hi on a straight segment, opens the 
system and the stadium's survival probability P{t) be- 
comes a useful statistical observable. Due to the bounc- 
ing ball orbits, P{t) is found to experience a cross-over 
from the above exponential decay (1) at short times to an 
algebraic decay ~ C/t at later times [24|. We point out 
here that the reason for this is that the stadium's clas- 
sical phase space is split into separate regions oc- 
cupied by 'fully-chaotic' and 'sticky' orbits, which 
are responsible for the exponential and algebraic 
decays respectively. As an orbit approaches the sticky 
region in phase space, which surrounds the bouncing ball 
orbits, it will inevitably escape quickly after it obtains 
an incidence angle \0\ < arctan (|^). This is a key point 
that will be discussed further in the two hole case shortly. 

We also remark that due to the splitting of the phase 
space, there is no justification for an intermediate purely 
exponential decay, as proposed generically for intermit- 
tent systems by Altmann et al. (see eq. (25) in [25[), but 
rather a coexistence of exponential and algebraic decay 
given by: 

T^f+\ _ /irregular, for t < i 
\e-^^ + f, fort>£, 

(3 In 3 + 4) ({a + /ij" )2 + (a - hf)^) 

n — V L CT] 

4(4a + 2^r)^ ' ^ ^ 

with parameters as defined in Figure JTj where we have 
neglected terms of order and i < can be found 
as described in |2^ and guarantees the splitting of the 
phase space as described above. The 'irregular' short- 
time behaviour is a result of geometry dependent short 
orbits which become less important if the hole is small. 
We note that the coefficient of the exponential term in (2) 
is 1 since for small holes and times greater than ^ 1/A, 
mixing causes the system to forget its initial state and 
therefore the probability decays as a Poisson process. 

Consider now the case of the stadium with two holes 
as shown in Figure [T] In Figure |2]we plot in the top 
panel a picture of the phase space, showing in different 
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FIG. 2: (Colour online) Phase space of open stadium with 2 
holes. Top: Initial conditions which will escape through hole Hi 
are shown in light yellow while those escaping through hole H2 in 
dark blue. Bottom: Colour grading of initial conditions going from 
purple, to orange, to white corresponding to short, medium and 
long escape times, {a = 2^ r = 1^ hi = 0.5, = 0.25). 



colours, the different sets of initial conditions which even- 
tually exit through each hole. The bottom panel shows 
the time scales of escape as noted in the caption. We 
notice that the phase space is again separated, as de- 
scribed above, and that the sticky, long-surviving orbits 
escape only through the hole on the straight segment Hi . 
Restricting the initial density of particles to one of the 
holes defines the transport problem and establishes the 
schematic setup of quantum dots and microwave cavi- 
ties, where particles /waves are injected through one of 
the holes and allowed to escape through either, thus cre- 
ating a direct link with experiment. Looking at the spa- 
tial distribution of the final (escape) coordinates (x/,6>/) 
(see Figure [s]) we observe that the long surviving orbits 
entering and subsequently exiting through i^i, unlike in 
the other possible entry/exit combinations, accumulate 
on the edges of the hole Xf = hf ^ 5^ ((5 <C 1) and have 
small angles Of. Note that (xf^Of) (hf^O^) as the 
time of escape tf ^ oo. This further confirms the split- 
ting of the phase space, but also that the classical spatial 
distribution of exiting particles has a well defined time- 
dependent character, which only exists in the situation 
described and plotted in Figure [Sj 

We define transmission and refiection survival proba- 
bilities by P- (t) and P-{t) respectively (z, j = 1,2), such 
that 

P^{t) = P{xi .,.xx^H\xoeHi,XfeHj), (4) 

where H = HiU M{xQ^t) is the number of collisions 
with the boundary up to time t and Xn denotes the po- 
sition of the particle at the nth collision. Hence, Pi if) 
is the probability that a particle injected from hole Hi 
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FIG. 3: (Colour online) 3D plot of the final (escape) coordinates 
and time of escape {xf,6f,tf) for the case of entry and exit through 
Hi. Only in this case are the 2 dark spikes observed. For other 
entry/exit combinations, a uniform 'carpet' with an exponentially 
decaying number of particles for larger t would be observed. The 
colour grading emphasizes the magnitude of the survival time of 
each orbit, {a = 2, r = 1, hi = 0.2, /i+ = 0.1). Ttaii ^ 631.85 is 
explained in Figure ^ 



will survive until time t given that it will escape through 
hole H2. It follows that only Pl(t)^ out of the four pos- 
sible survival distributions, has an algebraic tail, while 
the other three decay exponentially with an escape rate 
given by 7 = ^y^- Also, for Pi (t), at t > 1 and 6> <C 1, 
only particles starting near the edges of the hole are not 
immediately reflected back into the hole. Therefore, the 

extra constraint \6\ > arctan (± depending on 

the sign of 0) gives Pi{t) an algebraic tail 0{t~'^)^ as 
expected in integrable scattering problems. In summary 
the total survival probability Pi(t), where the subscript 



i indicates the injecting hole, is given by: 



^i(^) 



-ft 



D 



Pity 



-It 



^^^ple-^\ (5) 



pi{t) 



pl{t) 



(6) 



for t > t, where the are the respective asymptotic 
{t 00) reflection and transmission coefficients. Notice 



that p\ 



p1 



1 due to flux conservation, and pf = pi 



due to time-reversal symmetry. D is given by a similar 
calculation to [26^ : 

r(3 In 3 + 4) ({a + /ij" )2 + (a - hf)^) 
D= ^ ^ ^. (7) 

In Figure [4j we plot the four conditional distribu- 
tions P- (t) as functions of time t, and find an excel- 
lent agreement with the analytical results summarised in 
equations (5-6). We emphasize that the power law decay 
of Pi{t) is due to the geometric asymmetry of the hole's 
positions which exploit the marginally unstable bouncing 
ball orbits as to force a preference of escape through Hi . 
Furthermore, the splitting of the phase space into fully 
chaotic and sticky regions renders the later inaccessible 
to particles injected through H2. This would not have 
been be the case if both Hi and H2 were placed on a 



straight or curved segment. This important observa- 
tion offers the simplest possible example where 
a classically fully chaotic system exhibits time- 
asymmetric transport when opened. This is ex- 
pected to be relevant to many other intermittent systems. 
Also, the variety of options with regards to hole positions 
and sizes and system parameters offers ways of calibrat- 
ing and controlling these classical distributions and hence 
encourages the possibility of experimental observation of 
the quantum analogue which we discuss next. 

At low temperatures (~ 15 mK), electronic trans- 
port through the gate electrodes (openings) of a 2D elec- 
tron gas (quantum dot) is ballistic [HI [12]. For typi- 
cal semiconductor nano-structure parameters, the time 
scale Ttaii at which the above observed algebraic tail be- 
comes visible (see Figures [s] and [4| is of the order of a 
nanosecond (assuming an electron speed v ^ 10^ ms~^). 
This is slightly larger than the predicted Ehrenfest time 
te = \n{l/h) for chaotic systems [8 (the time scale 
at which quantum interference effects become apparent 
^ 0.3 ns), and thus at first instance suggests that direct 
observation of a quantum difference in transmission and 
reflection survival probabilities is unlikely in existing de- 
vices. We find that by varying the size and hole positions 
of the dot (while remaining in the ballistic regime) it is 
possible to calibrate and reduce Tfaii by a whole order of 
magnitude. A good way to do this is by elongating the 
stadium slightly such that a/r ^ 5 and by placing Hi 
at the very edge of the straight segment. However, since 
the nature of chaos lies in orbital instability, the Ehren- 
fest time Te calculated from the average Lyapunov 
exponent A does not faithfully represent quantum 
spreading of the near-bouncing ball orbits, which 
are responsible for the algebraic decay of Pl{t) noted 
here. For these, the finite-time local Lyapunov expo- 
nent is zero [27], therefore leading to a much longer va- 
lidity and persistence of the classical description of the 
sticky region in phase space. In fact, this region could be 
thought of as an hi depended fictitious island of stabil- 
ity in which loss of quantum-to-classical correspondence 
is much slower, resembling that in mixed systems, such 
that Te oc h~^^^ [28]. Therefore, experimental realiza- 
tion of the asymmetric transport (5) and (6) would not 
only confirm classical-quantum correspondence of 
the time-scales proposed in this particular sys- 
tem, but also emphasize the need for quantum 
corrections for transport in intermittent chaotic 
systems. 

Suppose we apply a time-dependent voltage V{t) 
across the gates of the stadium heterostructure such that 
the incoming current /*^(t) through hole Hi is propor- 
tional to V{t). Then the charge exiting through each hole 
will follow the driving current with a lag-time r which is 
distributed according to (5) or (6) appropriately. This 
can be modeled by 



dPUr) 



(8) 



where i and j indicate the injecting and exiting hole re- 
spectively. The observed, net current through the system 
is thus given by /^^^(t) = /r(^) + ^i(^) + ^^{t). Because 
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FIG. 4: Slightly offset plots comparing numerical simulations of 
the conditional survival probabilities (empty circles and squares) 
with the analytic expressions (5) and (6) as functions of time t in ns. 
The simulations consist of 10^ particles with stadium parameters 
given by: a = 2 jim^ r = 1 jim^ hi = 0.2 /^m and = 0. rtaii ~ 
6.315 ns is the large solution of e~^* = P2 7 where p\ ~ 0.5594 
was calculated numerically. 



the probability density ^^^^^^ is slightly skewed to the 
right, relative to the other densities, the two observables 
/f^*(t) and I^^\t) win differ by 

Plf {Pl{t)-Pi{t))dr. (9) 

For experimental observation we propose using a square 
wave signal V(t) = (1 + sign{smujt)) such that u > 
^/^taii as to accentuate the power-law contribution of 
Pi{t). Quantum interference effects such as universal 
conductance fluctuations may be statistically removed 

since the skewness of ^^J^^ is to leading order geometry 
dependent through the constant D in (7). In experiments 



of course, one should make sure that the excess density 
of charged particles within the dot is always low enough 
as to avoid a build up of an internal electric field which 
would effectively destroy the fictitious island of stability 
(sticky region) enclosing the bouncing ball orbits. For 
microwave billiards this is not an issue. 

To conclude, we have investigated the classical dynam- 
ics of the chaotic stadium billiard with one and two holes, 
and have found that the transmission and reflection sur- 
vival distributions of the latter case can have algebraic 
and exponential decays observed in the same classically 
ergodic geometry. We have identified the reason for this 
being the hole's asymmetric positioning on the straight 
segment of the billiard, which essentially splits the clas- 
sical phase space of the system, rendering the sticky re- 
gion surrounding the bouncing ball orbits inaccessible to 
chaotic orbits. As a result, the transmission and reflec- 
tion survival distributions are qualitatively different. We 
have obtained analytic expressions, confirmed them nu- 
merically and propose that observation of this classical 
result in semiconductor nano-structures (quantum dots) 
or microwave cavities would improve our understand- 
ing of classical to quantum correspondence in transport 
problems. Specifically it would imply that the Ehren- 
fest time and more generically quantum chaos predictions 
have correction terms subject to the underlying classical 
dynamics of the corresponding open systems. Further- 
more, this would possibly introduce new ways of calibrat- 
ing and controlling transport through complex systems 
by utilizing the sticky (non-mixing) channels of intermit- 
tent systems. 
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